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1. Introduction 

Any continuous physical model is empirically equivalent to a certain finite model. 
This thesis is widely used in practice: solution of differential equations by the finite 
difference method or by using truncated series is typical example. 

It is often believed that continuous models are “more fundamental” than discrete 
or finite ones. However, there are many indications that nature is fundamentally 
discrete at small (Planck) scales, and is possibly finite. 1 Moreover, description of 
physical systems by, e.g., differential equations can not be fundamental in principle, 
since it is based on approximations of the form / (x) ~ / ( x 0 ) + V/ ( x 0 )Ax . 

This essay advocates the view that finite models provide a more relevant de¬ 
scription of physical reality than continuous models which are only approximations 
in the limit of large numbers. 2 Using simple combinatorial models, we show how 
such concepts as continuous symmetries, the principle of least action, Lagrangians, 
deterministic evolution equations, etc. arise from combinatorial structures as a re¬ 
sult of the large number approximation. We also consider some approaches to the 
construction of discrete models of quantum behavior and related models describing 
evolution of gauge connections. 

Any statistical description assumes one or another concept of macrostate. We 
define macrostates as equivalence classes of microstates. This definition is especially 
convenient for models incorporating symmetry groups. We distinguish two types 
of statistical models: 

1. Isolated system is purely combinatorial object in the sense that “probability” 
of a microstate has a priori nature. Namely, all microstates are equiprobable, 
so their probabilities are equal to the inverse of their total number. The 
macrostates are specified by an equivalence relation on microstates. 

2. Open system is obtained from an isolated system by the following modifi¬ 
cation: The macrostates are specified by the same equivalence relation, but 
the probabilities of microstates depend on some parameters, which are in¬ 
troduced for approximate description of interaction of the system with the 
environment. 

The archetypal examples of isolated and open systems are, respectively, microcanon- 
ical (macrostates are defined as collections of microstates with equal energies) and 
canonical (macrostates are defined similarly, and interaction with the environment 
is parameterized by the temperature) ensembles. 

The classical description of a (reversible) dynamical system looks schematically 
as follows. There are a set X of states and a group G c \ < Sym (X) of transfor¬ 
mations (bijections) of X. Evolutions of X are described by sequences of group 

1 The total number of binary degrees of freedom in the Universe is about 10 122 as estimated 
via the holographic principle and the Bekenstein Hawking formula. 

2 Comparing the Planck length, ~ 10 -35 m, with the minimum length observable in experi¬ 
ment, ~ 10~ 15 m, we may assume that the emergence of the empirically perceived continuous 
space is provided by averaging over about 10 20 discrete elements. 
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elements g t G G c i parameterized by the continuous time t G T = [t a ,td\ C M. The 
observables are functions h : X —)■ M. 

We can “quantize” an arbitrary set X by assigning numbers from a number 
system T to the elements x G X, i.e., by interpreting X as a basis of a module. 
The quantum description of a dynamical system assumes that the module associated 
with the set of classical states X is a Hilbert space Hx over the held of complex 
numbers, i.e., T = C; the transformations g t and observables h are replaced by 
unitary U t and Hermitian H operators on H\ . respectively. 

To make the quantum description constructive, we propose the following mod¬ 
ifications. We assume that the set X is finite. Operators U t belong to the group 
of unitary transformations of the Hilbert space Aut ('Hx). Using the fact that this 
group contains a finitely generated — and hence residually finite — dense subgroup, 
we can replace Aut (Hx) by unitary representation of some finite group G that is 
suitable to provide an empirically equivalent description of a particular problem. 
A plausible assumption about the nature of quantum amplitudes implies that the 
held C can be replaced by an abelian number held 1 T . This held is a subheld 
of a certain cyclotomic held Q m which in turn is a subheld of the complex held: 
T < Q m < C. The natural number m, called conductor , is determined by the struc¬ 
ture of the group G. Note that the helds Q m and C provide empirically equivalent 
descriptions in any applications, because Q m is a dense subheld of C for any m > 3. 

In this paper we will assume that the time T is discrete and can be represented 
as a sequence of integers, typically T = [0, 1 ,..., T], 

Note also that the subscript a in the notation H a for Hilbert spaces is overloaded 
and can mean, depending on context: dimension of the space, a set on which the 
space is spanned, a group whose representation space is H a , etc. 

2. Scheme of statistical description 

For convenience of presentation, let us hx some notation: 

• U is the full set of states of a system A. The states from U are usually called 
“microstates” in statistical mechanics. 

• N = \U\ is the total number of microstates. 

• p u is the probability ( weight ) of a microstate u G U, p u = 1. 

u&U 

• ~ is an equivalence relation on the set U. 

• Taking an element u G U, we dehne the macrostate A as the equivalence class 
A = {v G U | v ~ u}. 

• A denotes the set of macrostates. 

: An abelian number field is an algebraic extension of Q with abelian Galois group. According 
to the Kronecker-Weber theorem, any such extension is contained in some cyclotomic field. 
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• The equivalence relation ~ determines the partition U — ]J A. 

aga 

• K — |A| is the number of macrostates. 

• N\ = | A| is the size of a macrostate A G A. 

• P\ denotes the probability of an arbitrary microstate from U to belong to the 
macrostate A. 


Isolated systems. The probability of a microstate of an isolated system is defined 
naturally 1 as p u = 1/N for any u G U, and the probability of any microstate 
from U to belong to the macrostate A is, respectively, P\ = N\/N. Since the 
“probabilities” in isolated systems have an a priori nature, such systems are in fact 
purely combinatorial objects, and we can talk about the number of combinations 
instead of probability. 


Open systems interact with the environment. This interaction is parameterized 
by assigning, in accordance with some rule, probabilities to all microstates. That 
is, the probability of a microstate uG A is a function 

Pu = Pu(o' 1 ,a 2 ,. ..) ( 1 ) 

of some parameters ai,« 2 ,.. .. These parameters and function (1) are determined 
by the specifics of a particular problem. For open systems P\ = ’Yh u &\Pu ( a i> a 2 , • • •)• 


Entropy. One of the central issues of the statistical description is the search for 
the most probable macrostates, i.e. the macrostates with the maximum value of 
P\. Technically, entropy is defined as the logarithm of the number (or probability) 
of microstates that belong to a particular macrostate. The concept of entropy is 
convenient for two reasons: 

1. Since the logarithm is a monotonic function, the logarithm of any function 
has the same extrema as the function itself. 

2. If a system can be represented as a combination of two independent systems, 
A = A'® A", then the macrostates of A can be represented as A'® A". So, when 
computing entropy of such decompositions, we can replace multiplication by 
a simpler operation — addition: log (PyPy) = logPy + logPA". 


1 This is “the equal a priori probability postulate” of statistical mechanics [1]. 
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Stirling’s formula is one of the main tools for obtaining continuum approxima¬ 
tions of combinatorial expressions: 


Inn! ~ 


n In n — n + - In (27rn) + 

superlinear s * v ^ 
logarithmic 


1 1 

12n ~~ 360n 3 + " ' 

V -V-' 

decreasing terms 


( 2 ) 


For our purposes it is sufficient to retain only the terms which grow with n, i.e. the 
superlinear and logarithmic terms. 


2.1 Examples of isolated and open systems 

To illustrate the above, let us give a few examples of isolated and open systems. 


Sequences of symbols. Isolated system. Let £ be an alphabet of size M: 

(oy,..., <J m , • • •, o"m) • (3) 

The microstates are sequences of the length T of symbols from £: 


u = a\ ■ ■ ■ a t ■ ■ ■ ar € U, at G £. 

The total number of microstates is N — M T . Let k u = (fc“,..., k^, ..., k^) be 
the vector of multiplicities of symbols a m in a microstate u. It is obvious that 
&“ + •■• + k 1 ^ = T. We define the equivalence relation ~ as follows 

u ~ v k u = k v = k. (4) 


The macrostate A defined by equivalence (4), consists of all sequences with the 
multiplicity vector k. The total number of macrostates is 

T + M-l 
M-l 

The size of the macrostate is 

T\ 

Afe ki\k 2 \ ■ ■ -k M \ 

Introducing the vector of “frequencies” (/i = k\/T,... Dm = k^i/T) and apply¬ 
ing the leading part of Stirling approximation log n! ~ n log n — n to the entropy 
S\ k = logN Xk of the macrostate A*, we obtain S\ k ~ TH(X), where 

M 

i= 1 

is the Shannon entropy [2] of a random variable X whose M outcomes have prob¬ 
abilities fi,, f M . 
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The model of symmetric random walk [3] is a slight modification of the 
above isolated system. Alphabet (3) contains now an even number of symbols 
M = 2 d, and it is divided into two parts £ = £ + ]j £_, where 

£ + = {(Tj G £ | 1 < i < d} and £_ = {a* G£|d+l<i< M} . 

The elements of £ + and £_ can be interpreted, respectively, as “positive” and 
“negative” unit steps in the d directions of the integer lattice Z d . 

The vector of multiplicities of symbols from £ can be written as k = (k + , k_), 
where k + = (k±,..., kd) and = (kd+ 1 ,..., k 2 d)- The equivalence of microstates u 
and v is defined now as follows 


u 


h u — h v — h v 

Ay - Ay | A/ 


( 5 ) 


The partition defined by (5) is a coarsening 1 of the partition defined by (4). The 
numbers of equivalence classes are figurate numbers of some d-dimensional regular 
convex polytopes — d-dimensional analogues of the octahedron. For example, in 
the case d = 3 the number of macrostates is equal to the (T + l)th octahedral 
number: 


K = (T + 1 ) 


2 (T + l) 2 + 1 
3 


Multinomial distribution. Open system. The microstates are also sequences 
of symbols from alphabet (3). But now the presence of an environment is assumed. 
The influence of the environment is parameterized by the assumption that any 

M 

symbol a m G £ comes with a fixed individual probability a m , such that a m = 1- 

m= 1 

Thus, oq 1 ^ 2 • • • = p\ k is the probability of a microstate from a macrostate A*,, 

defined by equivalence relation (4). The probability of a microstate from U to 
belong to the macrostate is described by the multinomial distribution: 


Microcanonical ensemble. Isolated system. The concept of a microcanonical 
ensemble is based on the classification of microstates by energy. More specifically, 
if there is a real-valued function on microstates E :[/—>■ M, then we can impose 
the equivalence relation on U: 

u ~ v E (u) — E (v) = E, u,v G U. (7) 

1 The partition V 2 = { Bi} of a set S is a coarsening of the partition V\ = {A,} of the same set 

if for every subset Aj G Vi there is a subset Bi £ V 2 such that Aj C Bi. The opposite relation 
among partitions is called the refinement [2]. 
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Assuming that the energy E takes a finite number of values: E e {E 1; E 2 ,..., E^} , 
we define the microcanonical ensemble as the macrostate A& which is an equivalence 
class of relation (7), i.e., the set of microstates with the energy E*,. More formally: 


Afc = {u | u G U A E [u) = E fc } . 


( 8 ) 


In statistical mechanics [4] the microcanonical ensemble is defined in terms of the 
Boltzmann entropy formula 

S Xk = k B In N\ k 

or, equivalently, via the microcanonical partition function 

N Xk = e 5 V kfl , 


where N Xk = | A* | and is the Boltzmann constant (we may assume that k B — 1). 


Canonical ensemble. Open system. A canonical ensemble is an open coun¬ 
terpart of the microcanonical ensemble. Macrostates of the canonical ensemble are 
defined by (8). Probability (1), that parameterizes the interaction with the envi¬ 
ronment in our general scheme, is now a function of a single parameter a = T, the 
temperature of the environment. Namely, the probability of a microstate u E is 
given by the Gibbs formula 


Pa. (T) 



E*/k „T 

5 


where the normalization constant 


z =iz N Xk e ~ Ek/kBT 

k =1 

is called the canonical partition function. 

3. Continuum approximations 

In this section we show that concepts such as continuous symmetry and the prin¬ 
ciple of least action may be obtained from combinatorial models as a result of the 
transition to the limit of large numbers. As an illustration, consider the open sys¬ 
tem described by distribution (6). In the case M = 2 ( binomial distribution), all 
calculations can be done explicitly. In this case both equivalence relations (4) and 
(5) coincide in virtue of the equality k\ + = T. Entropy of (6) for M = 2 takes 

the form 


S = In T\ — In kd — In {T — ki)\ + k\ lnaq + (T — kf) In a 2 . 


( 9 ) 
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Applying the growing terms of Stirling’s approximation (2) to this formula we have 

S ~ ^Stirling ‘S'superlin ~h ^log; where 


‘S'superlin T\nT k\ 111 



(T — ki) In 


T-k i 
«2 


( 10 ) 


and 

If T 

los 2 \2'Kki{T — ki) 

are superlinear and logarithmic parts of Stirling, respectively. Since S\ og and its 
derivatives are small for large T and ki, the maximum of entropy (9) is close to 
that of its leading part (10). Thus, the approximate extremum point fc* = a±T is 
obtained by solving the equation 


dS. 


superlin 


dk\ 



+ In 


T — k\ 
«2 


k\ 

T-k i 


cq 

1 — aq 


Further, we can expand Stirling around the point k\. Retaining terms up to the 
second order, we obtain the chain of approximations 


S ~ S Stirling ~ S* 


S, 


ld 2 S 


log I ki=a\T 


+ 77- 


superlin 


dk\ 


x (Jbi — a±T) 2 , 

k\=Oi\T 


which leads to the final formula 


S^S* = In 


27r7 1 CK 1 a2 2 Ta 1 a 2 


(h - a x Tf 


( 11 ) 


3.1 On the origin of continuous symmetries 

It is well known [3] that the large numbers asymptotic of the probability distribution 
of symmetric random walk on the lattice 7L d is the fundamental solution of the heat 
equation, also called the heat kernel 


K (t, x ) 


1 

(4tt t) d/2 


exp 


f x\ + xl + --- + x 2 d \ 

V 4 1 )■ 


( 12 ) 


Here t G M+ and Xi G M are the continual substitutes for T G N and for the 
difference ki — kd+i G Z, respectively. 

This gives an example of the emergence of continuous symmetries from the 
large numbers approximation. The symmetry group of the integer lattice Z d has 
the structure of the semidirect product 7L d x Gd- For simplicity, we can drop the 
normal subgroup Z d , the subgroup of translations, as inessential for our purposes. 
The group Gd is isomorphic to the semidirect product (Z 2 ) d x or, equivalently, 
to the wreath product Z 2 l S d- The size of Gd is equal to 2 d d\. For example, for 
the square lattice the group Gd =2 is the symmetry group of a square — dihedral 
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group of order 8. On the other hand, approximate expression (12) is symmetric 
with respect to the orthogonal group 0(d, R) with cardinality of continuum. 

The Lorentz symmetries — at least in 1+1 dimensions — can also be obtained 
in a similar way. Let us consider approximation (11) for the entropy of binomial 
distribution. We introduce the following continuous substitutes 


x := ki — k 2 , 
t:=T, 

v := — a 2 - (13) 


Obviously, — 1 < v < 1. With these substitutions, the approximation of binomial 
distribution takes the form 


P* (x, t ) = e s * 


7T (1 — V 2 ) t 


exp 



x — vt \ 
ffl — v 2 ) 



(14) 


The continuous variables x, t, and v may be called, respectively, the “space”, “time”, 
and “velocity”. 1 Expression (14) is the fundamental solution of the equation 

dP* (x,t) dP* (x,t) (1 — v 2 ) d 2 P* (x, t) 

dt V dx 2 dx 2 

This equation is called — depending on interpretation of the function P* (x, t ) - 
the heat, or diffusion, or Fokker-Plank equation. In the “limit of the speed of light” 
|w| = 1 equation (15) turns into the wave equation 


dP* (x, t ) dP* (x, t ) 

dt dx 


Let us introduce the change of variables: t = T H + t' and x = vT H + x' . If we assume 
that t' <C T h (i.e., T H can be thought as a “Hubble time”, and t! as a “typical time 
of observation”), then (14) can be rewritten as 


P* (x,t) 


/ 


2 

tt(1 -v 2 )T h 


exp 


1 ( x> — vt' 

2 T h V ffl — v 2 



The principal part of this expression is “relativistically invariant”. 


3.2 The least action principle as the principle of selection of the 
most likely configurations 

Let us compare the exact probability distributions with their continuum approxi¬ 
mations within individual equivalence classes of relation (5). 

1 In the paper [5], which is devoted to the “Zitterbewegung” effect in the 1 + 1 dimensional 

Dirac equation, a “drift velocity” is defined — just like in (13) — as the difference of probabilities 
of steps in opposite directions. It is shown that this definition leads to the relativistic velocity 
addition rule: w = (u + v) / (1 + uv). 
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Exact distributions. In the binomial case, an equivalence class of (5) is defined 
by fixing the difference k\ — k 2 X. We denote the equivalence class of sequences 
connecting the space-time points (0, 0) and (X, T) by A x,t- The size of A x,r is 
equal to 

T\ T\ 

Nx - t = ww " mwy ■ 

The binomial distribution in terms of the variables A" and T and parameter v takes 
the form 


P(X,T) = 


T\ 


m !(^)! 


1 + V 


T+X 

2 


Consider an increasing sequence of time instants (“times of observations”) 


r — {T 0 — 0,..., Tj_ 1; Tj,..., T n — T} , T t -\ < Tj. 


(16) 


Let us select trajectories that pass through the sequence of spatial points 


X — {X 0 — 0,..., X n — A} 


(17) 


corresponding to sequence of times (16). Admissible trajectories must satisfy the 
inequality \Xi — X,_i < Tj — T i _ 1 — “the light cone restriction”. According to the 
conditional probability rule, the probability of the trajectory (y, r) is equal to 


p = 

1 X,T 


P (X, T) p 


iPW-Xt-uTf-Ti-t) 


i =1 


(¥)!(¥)! 


T\ 


n 

i— 1 


(Ti — Tj-i)! 


Tj+Xj _ Tj-i+Xj-! \ , ( T t -X x _ Tj—i—Xj-i 
2 2 / V 2 2 


For a given sequence of time instants (16) one can formulate the problem of find¬ 
ing trajectories with maximum probability P T = max P x T . This can be done by 

x 

searching among all admissible sequences (17). In the general case, there are many 
distinct trajectories with the same maximum probability, i.e., we do not have here 
a “deterministic” trajectory. 


Continuum approximation. To apply our reasoning to approximate distribu¬ 
tion (14), we will consider the time sequence 


T "(tg t a , . . . , tj—1, tj, • • • , C 


4 } 


together with respective sequences of spacial points 


X = {x 0 = x a , ..., Xi- 1 , Xi,...,x n = x b } . (18) 

We assume that the time points are equidistant: ti — ti -1 = At, and we will use the 
notation Aaq = x* — aq_ i. 
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Now the approximate probability of a trajectory (y, r) that connects the space-time 
points (x a ,t a ) and (x b ,t b ) takes the form 


p* _ 


exp 


1 / X b - X a ~ V (t b ~ t a ) 


7T (1 - V 2 ) (t b -t a ) 12 (t b ~ t a ) V VT 3 


X 


Qx 

A n 


where 


and 


A = 


7T (1 — V 2 ) At 


n 

Qx,r = n exp 

i=1 


1 / A Xi — vAt\ 2 1 

2A t V J J 


= exp 



Axi/At — v\ 

y/1 — V 2 J 



(19) 

( 20 ) 


The summation of factors in (19) over all values Ax* G (—oo, oo) reproduces correct 
normalization of probabilities for any time slice: 


OO 



exp 


—oo 


1 / Ax* — uAt\ 2 | d(Axi) 

2A t V ) J 2“ 


/ 


7T (1 — V 2 ) At 
2 


= A. 


Note that this normalization 1 is an approximation which is incompatible with the 
“speed of light limitation”: — At < Ax^ < At. 

Replacing the sequence of spacial points (18) by a differentiable function x (t) 
such that x (t a ) = x a , x ( t b ) = x b , introducing approximation A Xi ~ x (t) At and 
taking the limit n —y oo we can write instead of (20) the formula 

Qx,r ~exp|-^[x(f)]|, 


where 

S[x(t)] = J Ldt = J ($=£) it. 

ta ta 

Assuming for a while that v depends on t and x, we obtain the following Euler- 
Lagrange equation 


d dL dL 
dt dx dx 


, 9 x o dv dv dv 

0 =>■ x(l — v)+xv— — [~2xv— v — 

v ' ox at ox 


(i+» 2 ) 


dv 

dt 


0. 


X A similar normalization is one of the cornerstones of the path integral formalism [6]. 
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Clearly, this equation describes “deterministic” trajectories. If we return to the 
initial assumption that v does not depend on the space-time variables, then the 
Euler-Lagrange equation reduces to the form 


x (t) = 0. 


This equation together with the boundary conditions gives the following formula 
for the extremals 

%b %a + ! %atb %bta 
' I * 

tb ta tb la 

i.e., the most probable trajectories are straight lines. 



4. Combinatorial models of quantum systems 

To build models that can reproduce quantum behavior, it is necessary to formulate 
the basic ingredients of quantum theory in a constructive way (for a more detailed 
consideration see [7]). 


4.1 Constructive core of quantum mechanics 

In traditional matrix formulation quantum evolutions are described by unitary op¬ 
erators in a Hilbert space "H. Evolution operators U belong to a unitary repre¬ 
sentation of the continuous group Aut fiH) of automorphisms of "H. To make the 
problem constructive we should replace the group Aut fiH) by some finite group G 
which should be empirically equivalent to (a subgroup of) Aut fiH). 

The theory of quantum computing [8] proves the existence of finite sets of univer¬ 
sal quantum gates that can be combined into unitary matrices which approximate 
to arbitrary precision any unitary operator. In other words, there exists a finitely 
generated group which is a countable dense subgroup of the continuous group 
Aut (W). 

A group G is called residually finite [9], if for every g e G, g 1, there exists a 
homomorphism 0 from G onto a hnite group H , such that 0 {g) 1 . This means 

that any relation between the elements of G can be modeled by a relation between 
the elements of a hnite group. This can be illustrated by analogy with the widely 
used in physics technique, when an infinite space is replaced by, for example, a 
torus whose size is sufficient to hold the data related to a particular problem. 

According to the theorem of A.I. Mal’cev [10], every finitely generated group of 
matrices over any held is residually hnite. Thus we have the sequence of transitions 
from the group with cardinality of continuum through a countable group to a finite 


group: Aut (fi-i) 


approximation 


> G r 


homomorphism 


> G 


4.2 Permutations and natural quantum amplitudes 

As is well known, any linear representation of a hnite group is unitary. Any repre¬ 
sentation of a hnite group is a subrepresentation of some permutation representation 
(see, e.g., [11-15]). Let U be a representation of G in a K-dimensional Hilbert space 







Discrete dynamical models 


13 


Hk- Then U can be embedded into a permutation representation P of G in an 
N-dimensional Hilbert space Hn- where N > K. The representation P is equivalent 
to an action of G on a set of things = {uq,..., cjn} by permutations. In the 
proper case N > K, the embedding has the structure 



( 21 ) 


where 1 is the trivial one-dimensional representation, mandatory for any permuta¬ 
tion representation; V is a subrepresentation, which may be missing. T is a matrix 
of transition from the basis of the representation P to the basis in which the permu¬ 
tation space "Him is split into the invariant subspaces Hn-k an d Hk- Evolutions in 
the spaces Hk and Hn-k are independent since both spaces are invariant subspaces 
of Hn- So we can treat the data in Hn-k as “hidden parameters” with respect to 
the data in Hk- 

A trivial approach would be to set arbitrary (e.g., zero) data in the complemen¬ 
tary subspace Hn-k- This approach is not interesting since it is not falsifiable by 
means of standard quantum mechanics. In fact, it leads to standard quantum me¬ 
chanics modulo the empirically unobservable distinction between the “finite” and 
the “infinite”. The only difference is technical: we can replace the linear algebra in 
the K-dimensional space Hk by permutations of N things. 

A more promising approach requires some changes in the concept of quantum 
amplitudes. We assume [7, 16,17] that quantum amplitudes are projections onto 
invariant subspaces of vectors of multiplicities of elements of the set on which 
the group G acts by permutations. The vectors of multiplicities 



n) 


( 22 ) 


\"N/ 


are elements of the module Hn = N N , where N = {0,1,2,...} is the semiring of 
natural numbers. Initially we deal with the natural permutation representation of G 
in the module Hn- Using the fact that all eigenvalues of any linear representation of 
a finite group are roots of unity, we can turn the module H N into a Hilbert space Hn- 
It is sufficient to add Cth roots of unity to the natural numbers to form a semiring , 
which we denote by Nc- The natural number C, called conductor, is (a divisor of) 
the exponent of G, which is defined as the least common multiple of the orders of 
elements of G. In the case C > 2 the negative integers can be introduced and the 
semiring Nc becomes a ring of cyclotomic integers. To complete the conversion of 
the module H N into the Hilbert space Hn, we introduce the cyclotomic field Q c as 
a field of fractions of the ring Me- 1 If C > 3, then Qc is a dense subfield of the field 

: By taking into account symmetries of a specific problem, we can use instead of Qc some its 
subfield, an abelian number field T < Qc- 






14 


V. V. Kornyak 


of complex numbers C. In fact, algebraic properties of elements of Qc are quite 
sufficient for all our purposes — for example, complex conjugation corresponds to 
the transformation r k — y r c ~ k for roots of unity — so we can forget the possibility 
to embed Qc into C (as well as the very existence of the field C). 

Thus, we will assume that "Hk in decomposition (21) is a Hilbert space over the 
field Qc, and quantum amplitudes are elements of "Hk of the form \fj) = Pru |n), 
where Pru is the projection operator from the module Hn onto the space "Hk cor¬ 
responding to the subrepresentation U in (21). 

4.3 Measurements and the Born rule 

The general scheme of measurements 1 in quantum mechanics is reduced to the 
following. 

• A partition of the Hilbert space into mutually orthogonal subspaces is given: 

U = Ui@ ■■•©Hi®-- - . 

Typically jii are eigenspaces of some Hermitian operator A (i.e. A = Af) 
called an “observable”. 

• There is a measuring device configured to select a state \<f) G "H, of a quantum 
system. 

I Yes, measuring device responds; 

• A result of a single measurement = < 

I No, no response. 

In accordance with the projection postulate , the output Yes is interpreted as 
transition of the system into the state \<f>) after the measurement. 

If jii is an eigenspace of an observable A with eigenvalue a, it is said that the 
“outcome of the measurement is equal to” a. 

• Relative number of Yes in a set of measurements is described by the Born 
formula. 


The Born rule 2 states that the probability to register a particle described by the 
amplitude \if) by an apparatus configured to select the amplitude \f>) is 


P 


(<m)l 2 

{<t> ' 


x To avoid inessential technical complications we consider here only the case of pure states. 

2 There have been many attempts to derive the Born rule from the other physical assumptions 
— the Schrodinger equation, many-worlds interpretation, etc. However, Gleason’s theorem [18] 
shows that the Born rule is a logical consequence of the very definition of a Hilbert space and has 
nothing to do with the laws of evolution of physical systems. 
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In the “finite” background the only reasonable interpretation of probability is the 
frequency interpretation: probability is the ratio of the number of “favorable” com¬ 
binations to the total number of combinations. So we expect that P(0, -0) must be 
a rational number if everything is arranged correctly. Thus, in our approach the 
usual non-constructive contraposition — complex numbers as intermediate values 
against real numbers as observable values — is replaced by the constructive one - 
irrationalities against nationals. From the constructive point of view, there is no 
fundamental difference between irrationalities and constructive complex numbers: 
both are elements of algebraic extensions. 

4.4 Illustration: natural amplitudes and invariant subspaces of per¬ 
mutation representation 

Consider the action of the alternating group A 5 on the vertices of the icosahedron. 
The group A 5 has a presentation of the form 

A 5 = (a, b | a 5 = b 2 = ( ab ) 3 = 1 ) . (23) 

The Cayley graph of this presentation is shown in Figure 1. 



Figure 1: Cayley graph of A 5 for presentation (23). Pentagons, hexagons and links 
between adjacent pentagons correspond to the relators a 5 , (ab) 3 and b 2 , respectively. 

A 5 has five irreducible representations: the trivial 1 and four faithful representations 
3,3', 4, 5; and three primitive 1 permutation representations having the following 
decompositions into the irreducible components: 5 = 1© 4, 6 = 1© 5, and 

10 ^ 1©4©5. 

: A transitive action of a group on a set is called primitive [13], if there is no non-trivial partition 
of the set, invariant under the action of the group. 
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The action of A 5 on the icosahedron vertices = {1,..., 12} is transitive, but 
imprimitive with the non-trivial partition into the following blocks 

(| B 1 | ■ ■ ■ | B t | ■ ■ ■ | B 6 |} = {| 1, 7 | ■ • ■ | i,i + 6 | ■ ■ ■ | 6, 12 |} , 

assuming the vertex numbering shown in Figure 2. Each block B\ consists of a pair 



9 


Figure 2: Icosahedron. Invariant blocks are pairs of opposite vertices. 

of opposite vertices of the icosahedron. Permutation representation of the action 
of A 5 on the icosahedron vertices has the following decomposition into irreducible 
components 

T~ x (12) T = 1 © 3 © 3' © 5, (24) 

where T is a matrix of transition from the “permutation” to the “splitting” basis. 

Actually there is no necessity to compute transformation matrices like T in 
(24) explicitly. There is a way [7] to express invariant scalar products in invariant 
subspaces in terms of easily computable matrices of orbitals [14,15], i.e., orbits of 
the action of a group G on the Cartesian product hi x hi. 

For the action of A 5 on the set of icosahedron vertices, the matrices of orbitals 
have the form 

- 4 ' = = ( 1 ° £)■•*= (y X ) ' = { Y x y) • < 25 > 

where I n is n x n identity matrix, and 


(0 

1 

1 

1 

1 

A 


/0 

0 

0 

0 

0 

°\ 

1 

0 

1 

0 

0 

1 


0 

0 

0 

1 

1 

0 

1 

1 

0 

1 

0 

0 


0 

0 

0 

0 

1 

1 

1 

0 

1 

0 

1 

0 

5 Y ~ 

0 

1 

0 

0 

0 

1 

1 

0 

0 

1 

0 

1 


0 

1 

1 

0 

0 

0 

V 1 

1 

0 

0 

1 

V 


\o 

0 

1 

1 

0 

0/ 
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In terms of matrices (25) the invariant bilinear forms (scalar products) correspond¬ 
ing to decomposition (24) take the form 


I3i — — (A-i + A.2 + -4.3 + .A 4 ), 


B3 — - [ - 4 i — - 4 2 


1 + 2r 2 + 2r 3 „ 1 + 2r 2 + 2r 3 „ 

--4s +--- A 4 


5 


5 


„ ~ . A 1 + 2r 2 + 2r 3 , 1 + 2r 2 + 2r 3 

13 3 ' — - ( Ai — A 2 H----43----44 

b,=^-(a 1 + a 2 -Ia 3 -^a 4 

12 V 5 5 


where r is a 5th primitive root of unity. It is easy to verify that the cyclotomic 
integer 1 + 2r 2 + 2r 3 is equal to — \/5. 

Let us consider the scalar products of projections of “natural” vectors. If 
projections of vectors with natural components m = (mi,..., mi 2 ) and n = 
(ni,..., ui 2 ) T onto the invariant subspaces corresponding to a = 1,3,3', 5 are 
and 'k Q , respectively, then (<L a | \k a ) = (m\B a \n). That is, we have 

($1 | vf 1 ) = -L (^(m\Ai\n) + (m\A 2 \n) + (m\A 3 \n) + (m\A 4 \n)^j , (26) 

($ 3 | ^ 3 ) = I ^(m|-4i|n) — (m|-4 2 |n) + -^ ^(m|«4 3 |ri) — (m|-4 4 |n))^ , 

($ 3 , | T 3 /) = i ((m\Ai\n)-(m\A2\n)-^^(m\A 3 \n) + (m\A4\n)^ , 

($5 | 4 5 ) = ^ r(m|A|n) + (m|-4 2 |n)-^(m|-4 3 |n) + (m|-4 4 |n)^ . 

Let us give two remarks on these expressions: 

• Scalar product (26) can be written as 

($1 | Ti) = -^(m 1 +m 2 ^ -f m 12 ) (ni + n 2 H- f n 12 ) • 

This is the general case: any permutation representation of any group con¬ 
tains the trivial one-dimensional subrepresentation with the scalar product 
like (26): 

<*> l *>> = * (|>) (|>)- 

The trivial subrepresentation can be interpreted as the “counter of particles”, 
since the linear permutation invariant n* is the total number of elements 
from 12 in the ensemble described by the vector n. 
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• The Born probabilities for subrepresentations 3 and 3' contain irrationalities 
that contradicts the frequency interpretation of probability for finite sets. 
Obviously, this is a consequence of the imprimitivity: one can not move an 
icosahedron vertex without simultaneous movement of its opposite. To resolve 
the contradiction, mutually conjugate subrepresentations 3 and 3' must be 
considered together. The scalar product 

($ 3 ffi 3 ' | ^ 303 ') = ^{m\Ai\n) - (m\A 2 \n)^ 

in the six-dimensional subrepresentation 3 © 3' always gives rational Born’s 
probabilities for vectors of multiplicities defined as in ( 22 ). 

4.5 Quantum evolution 

In standard quantum mechanics an elementary step of evolution is described by the 
Schrodinger equation 

I0> = H W) ■ 

In quantum mechanics based on a finite group G a step of evolution has the form 

\A+i) = u | ip t ), 

where U = U (g), g G G and U is an unitary representation of G. In this case, there 
is no need for a Hamiltonian, though, for comparison purposes, it can be introduced 
by the formula 

H — i In U = 7r(a 0 I + ot-JJ + ■ • • + a n _if/ n_1 ) , 

where II is the unit matrix; n is the period of U , i.e. U n = I ; aj 6 J are easily 
computable coefficients; 0 < k < n. The energy levels (eigenvalues) of H are 
2tt /c 

Ek = - . The non-algebraic (transcendental) number n appears here as the 

n 

result of summation of infinite series — the natural logarithm is essentially an 
infinite construct. 

Note that a single unitary evolution is physically trivial, as it describes only a 
change of coordinates (“rotation”) in a Hilbert space. Namely, for the evolution of 
a pair of vectors j 0 t) = U | 0 o), |V’t) = U l^o), we have 

(<MVr) = (0 0 \U ] U\ 0 O ) = (0o| V’o) • 

This means that a single deterministic evolution can not provide physically observ¬ 
able effects. Thus, a collection of different evolutions is needed. Suppose that the 
operators of evolution belong to a unitary representation U of a group G. Then 
two different evolutions U and V can be represented as U = U (g\g2 ■ ■ ■ gr) and 
V = U (/ 1/2 • • • / t ), where g±, ■ ■ ■ , gr', fi, ■ ■ ■ , /t G G. These evolutions provide a 
nontrivial physical effect if 

{<h I V’t) = <00 I V*u | 0o> ¥= <00 10o) 4 =^ V'U = U (/y 1 • • • f^ gi ■■■g T )^l 
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or 

h = It 1 - ••/fV- ■ -9t ~f~ 1, 

where 1 is the identity of G. The expression h = ffi 1 • • • f^gi ■ ■ ■ gr is called 
the holonomy at the point T of principal G-connection. In differential geometry, 
infinitesimal analogue of holonomy is called the curvature of the corresponding 
connection. As is well known, all fundamental physical forces are represented in 
the gauge theories as curvatures of appropriate connections. 

Observations support the view that fundamental indeterminism is really the 
modus operandi of nature. 1 In our view this indeterminism arises from a funda¬ 
mental impossibility to trace the identity of indistinguishable objects during their 
evolution. Hermann Weyl discussed this issue in detail in [19]. More formally, iden¬ 
tification of objects at different time points is provided by a connection (parallel 
transport). 

In our setting we consider dynamical system as a fiber bundle (W, T, fi, G, r) 
over discrete time T = [0,1,..., T], where typical fiber H = {cai,. ..,wn} is canon¬ 
ical set of states, structural group 

G = {gi,..., g m ,..., g M } < Sym (fi) = S N (27) 

is the group of symmetries of states, r is a projection W — > T. 

Connection ( parallel transport ) X ili2 e G defines isomorphism between the fibers 
at different times of observations: 

^<2 “ hi ' 

There is no objective way to choose the “correct” value for the connection. A 
priori , any element of G may serve as X tlt2 . In continuous gauge theories, the 
gauge fields (fields of connections) are determined from the principle of least action 
using Lagrangians chosen for different reasons. For example, in the case of the 
Yang-Mills theory (covering also the case of Maxwell equations), the Lagrangian 
Lym = Tr [F A *F] is used, where F is the curvature form of a gauge connection, 
* denotes the Hodge conjugation. Analysis of the structure of L yM in the dis¬ 
crete approximation [20] shows that it can be expressed in terms of traces of the 
fundamental representation of holonomies of a gauge group. 

4.6 Combinatorial models of gauge and quantum evolution 

Consider a simple combinatorial model involving random choice of the rules for 
identification of states of dynamical systems at different points of time. The time 
of our model is the sequence T = [0,1,..., t — 1, t, ..., T], A parallel transport 

1 There are persistent attempts to develop a deterministic version of quantum mechanics, as 
if determinism were a “synthetic a priori judgment” — an inevitable (though not deducible from 
logic alone) necessity. However, since the time of Kant up to now there are no convincing evidences 
of the very existence of judgments of this kind. More likely, the belief in determinism is a mental 
habit formed by long macroscopic experience. 
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connecting the initial and the final time points can be decomposed into the product 
of elementary steps: 

X 0 ,i • • • ''' X-t-i,t • (28) 

We assume that any elementary step is an element of group (27) with probability 
independent of time: 


M 

P (/^t—l,t g m ^ G) Oimi ^ ^ 1 • 

m= 1 

All possible paths (28) form the set of microstates U. The microstate corresponding 
to (28) is the sequence u = g mi ,.... g mt ,..., g mT , where g mt = X t _ u . Its probability 
is the product p u = a mi ■ ■ ■ a mt ■ ■ ■ a mT . 

We can dehne a natural equivalence relation on U as the triviality of the holon- 
omy of a pair of paths: 

A 0 ,i • • • X-t-i,t ~ Top • • • 7 t-i,t 7 r i ljT • • • 7 0 j j;X 0|1 • • • X T _ 1T = 1 . 

This equivalence allows us to define M macrostates Ai,..., Am- Statistical evolution 
of all the macrostates can be calculated simultaneously by a simple algorithm. The 
distribution of the macrostates at the moment T is the following element of the 
group algebra 


At = P\ 1 gl + P\ 2 & 2 + • • • + Ta m SM — («lgl + Oi 2g2 + • • • + CVMgM) • 


The algorithm of binary exponentiation computes this expression by performing 
O(logT) multiplications. In simple cases the probabilities P\ m can be written 
explicitly, e.g., for the cyclic group G = we have (assuming g x = 1) 


Px 


m 


E 

kl~\~k2~\ -l~^M—^ 

/C2“b2/C3H-h(M — 1)^M =772 


T\ 


k]\k 2 \ 


k\ k 2 
a 1 a 0 


■ ■ ■ a 


k M 

M ’ 


(mod M) 


Having a representation U of the group G in a Hilbert space T-L we can associate 
with gauge evolution (28) the quantum evolution: 

l^r) = u ''' X t-i,t ■ ■ ■ x oq) l^o), V’o, e U . 

Simulation of many important features of quantum behavior requires models 
involving spatial structures explicitly. The set of states of a system with space is 
the set of functions 

ft = E x , (29) 

where 

X = {xi,... ,x n } 
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is a space , and 

S = {(Ji, . . . , CT| S |} 

is a set of local states. Having the groups of spatial 

F={fi,...,f|F|}<Sym(X) 

and internal 

r = {7i) * • ■ > 7|r|} < Syrn(E) 

symmetries we can construct a symmetry group of the whole system . This group, 
having a structure of the wreath product 

G = T l x F = T x x F, 


acts on the set given by formula (29). 

It is worth to say a few words about the most common quantum models with 
spatial structures: quantum cellular automata [21] and quantum walks [22], It is 
proved that models of both types are able to perform any quantum computation, 
i.e., they can simulate quantum Turing machines. 

The Hilbert space of a quantum cellular automaton has the form 


nj _ o/0X 

tt — tt s , 


where X is usually a d- dimensional lattice: 7G l or its finite counterpart h d N (one can 
also take an arbitrary regular graph as a lattice X); "He is a Hilbert space associated 
with a set £ of local states of sites i6X. It is assumed that there is a local update 
rule U x , which is a unitary operator acting on the Hilbert space where JV X 

is a neighborhood of the point x. Since the neighborhoods of different points may 
intersect, some effort should be made to ensure global unitary. To provide the 
required compatibility several different definitions of quantum cellular automata 
were proposed. For properly defined automaton the local updates can be combined 
into a unitary operator U on "H that describes an elementary step of evolution of 
the whole system. Then the evolution of the system is defined by the operator U 7 . 

The spatial structure in a model of quantum walk is a /.-regular graph X. In 
the most usual case k = 2, the space X is taken to be either Z or 7L N . Let kix 
be the Hilbert space spanned by the vertices of X. The construction of a quantum 
walk uses also an auxiliary fc-dimensional Hilbert space "He, the “coin space”, and 
a fixed unitary “coin operator” C acting on ktc- A typical coin operator in the case 
k > 2 is the Grover coin ( Grover’s diffusion operator)-. 


G 


1 ^, where 



In the case of one-dimensional quantum walk (k = 2) many different one-qubit 
gates, like the Hadamard gate etc., are used. In particular, the Grover coin G 
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coincides with the Pauli-X gate at k = 2. The Hilbert space of the whole system is 
the product 

n = n x ® u c . 

Roughly speaking, the coin operator C “selects directions of spatial shifts”. The 
spatial shifts are performed by an unitary shift operator S acting on "Hx at con¬ 
ditions given by the coin C. T steps of evolution of the system are performed by 
the transformation U T , where U is the following combination of the coin and shift 
operators 

U = S (H fc <g> C). 


5. Summary 

Starting with the idea that any problem that has a meaningful empirical content 
can be formulated in constructive finite terms, we consider the possibility of deriva¬ 
tion of many important elements of physical theories in the framework of discrete 
combinatorial models. We show that such concepts as continuous symmetries, the 
principle of least action, Lagrangians, deterministic evolution equations can be 
obtained by applying the large number approximation to expressions for sizes of 
certain equivalence classes of combinatorial structures. 

We adhere to the view that quantum behavior can be explained by the funda¬ 
mental impossibility to trace identity of indistinguishable objects in the process of 
their evolution. Gauge connection is that structure which provides the identity: 
that is why the gauge holds are so important in quantum theory. 

Using general mathematical arguments we show that any quantum problem can 
be reduced to permutations. Quantum interferences are phenomena observed in 
invariant subspaces of permutation representations and expressed in terms of per¬ 
mutation invariants. In particular, this approach gives an immediate explanation 
for the appearance of complex numbers and unitarity in the formalism of quantum 
theory. 

We consider some approaches to the construction of discrete models of quantum 
behavior and related models describing evolution of gauge connections. 

Acknowledgements. I am grateful to A.Yu. Blinkov, V.P. Gerdt, A.M. Ishkhanyan 
and S.I. Vinitsky for many insightful discussions and comments. 
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